Stratified Graphene-Noble Metal Systems for Low-Loss Plasmonics Applications 
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We propose a composite layered structure for tunable, low-loss plasmon resonances, which con- 
sists of a noble-metal thin film coated in graphene and supported on a hexagonal boron nitride 
(hBN) substrate. We calculate electron energy loss spectra (EELS) for these structures, and nu- 
merically demonstrate that bulk plasmon losses in noble-metal films can be significantly reduced, 
and surface coupling enhanced, through the addition of a graphene coating and the wide-bandgap 
hBN substrate. Silver films with a trilayer graphene coating and hBN substrate demonstrated sur- 
face plasmon-dominant spectral profiles for metallic layers as thick as 34 nm. A continued- fraction 
expression for the effective dielectric function, based on a specular reflection model which includes 
boundary interactions, is used to systematically demonstrate plasmon peak tunability for a variety 
of configurations. Variations include substrate, plasmonic metal, and individual layer thickness for 
each material. Mesoscale calculation of EELS is performed with individual layer dielectric functions 
as input to the effective dielectric function calculation, from which the loss spectra are directly 
determined. 
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I. INTRODUCTION 

Plasmonic devices hold promise for a host of metama- 
terials applications, due to their ability to control light 
propagation on a subwavelength scale. Noble metal nano- 
materials are in some ways ideal components in such de- 
vices - they possess tunable, large amplitude plasmon 
resonances, which can be excited at optical wavelengths. 1 
Nevertheless, metallic plasmonic devices at infrared and 
visible wavelengths present significant challenges due to 
bulk plasmon losses. These large losses severely limit the 
practicality of these materials for a wide variety of ap- 
plications, particularly in telecommunications and pho- 
tovoltaicsi^ 

Bulk losses may be mitigated through the use of very 
thin noble-metal films; however, fabrication of uniform 
thin metal films is experimentally challenging. Thin films 
tend to form islands and often require adhesion layers 
that alter the electronic structure of the device 4^ Addi- 
tionally, thin noble-metal films still suffer significant re- 
sistive losses in the visible regime^ Composite materials 
designed to shift plasmon resonances to a low-loss regime, 
while increasing coupling to surface plasmons and dimin- 
ishing bulk resonances, provide clear advantages from a 
fabrication standpoint. Systems that combine the strong 
and tunable optical-wavelength plasmon resonances of 
noble metals with materials possessing improved trans- 
port properties could mitigate surface-plasmon losses in 
the visible regime. 

Graphene is an ideal candidate for such a composite 
structure due to its unparalleled carrier mobility. This 
allows for extremely enhanced and tunable electromag- 
netic response spectra when doped with other plasmonic 
materials, or fabricated as a component of a multilayer 
structure. Monolayer graphene has a response spectrum 
dominated by absorption peaks at « 4.5 eV and f=s 15 eV, 
the 7T and n + a surface plasmons. Graphene plasmon 
resonance provides low losses in the frequency regime be- 
low the optical phonon frequency of 0.2 eV, where large 



losses are typically present for metallic plasmonic mate- 
rials £ 

Hexagonal boron nitride (hBN) provides several advan- 
tages as a substrate for graphene-based plasmonic mate- 
rials, as it is an ultra-flat wide bandgap insulator with 
excellent thermal transport properties^ The flat mor- 
phology and uniformity provided by hBN are are desir- 
able as key factors in preserving the transport proper- 
ties of graphene.™ Results of our numerical study clearly 
elucidate the advantages of the hBN spectral response 
for low-loss plasmonics applications. We demonstrate 
in this work that through the combined advantages of 
graphene and noble-metal films, along with careful choice 
of substrate material and layer thickness, multilayer sys- 
tems can be tuned for low-loss surface plasmon resonance 
(SPR). 

Metamaterials derive their exotic properties in part 
from advantageous behaviors of their constituent ma- 
terials, as well as collective behaviors that emerge due 
to interactions within the system. The complex task of 
designing metamaterials tailored for applications as di- 
verse as photovoltaics, biosensing, and microscopy gives 
rise to a need to develop efficient methods for predicting 
metamaterial properties. These methods must realisti- 
cally treat both the individual material properties and 
the interactions among the constituent materials. We 
detail and employ such a method for the calculation of 
electron energy loss spectra (EELS) of multilayer struc- 
tures consisting of graphene layers on noble-metal (silver, 
gold, and copper) films with silicon and hBN substrates. 

The effective dielectric function is based on a specu- 
lar reflection model, first derived by Lambin et al£, and 
takes into account the boundary conditions across each 
layer in the stratified structure. The use of the efficient 
continued-fraction expression along with pre-prepared li- 
braries of dielectric functions for the individual materials 
opens up the possibility of multilayer graphene compos- 
ites by design. 
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B. The Effective Dielectric Function 



Graphene Layer(s) 



Noble Metal Layer 



Si Layer 
(Semi-Infinite) 
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Noble Metal Layer 



Vacuum 
(Semi-Infinite) 



Metallic Film 




Substrate: a) Si/Si0 2 b) hBN 



FIG. 1. Multilayer structure: EELS are calculated for a struc- 
ture consisting of graphene top layer(s), noble-metal middle 
layer, and semiconductor substrate. 



II. EELS CALCULATION DETAILS 
A. General Procedure 

EELS are calculated for a variety of multilayer sand- 
wich structures as depicted in Fig. [TJ In general, the 
sandwich structure consists of a noble-metal middle layer, 
a graphene coating of an integer number of layers, and a 
semiconducting substrate. Graphene layer numbers vary 
from 1 to 20. Noble metals employed include silver, gold, 
and copper, and are varied in thickness. The effect of 
two different substrates is considered: (1) a more tradi- 
tional substrate, Si with a 4 nm top layer of SiC>2 (Fig. Q] 
(a)), and (2) hBN (Fig. [1] (b)). 4nm is a typical thick- 
ness for the Si02 layer, which arises from the thermal 
processing of the Si substrate.— The Si layer is assumed 
to be semi-infinite. In the case of the hBN substrate, an 
integer number of hBN layers is considered, with a semi- 
infinite vacuum layer below the hBN layer. The single- 
layer thickness for both graphene and hBN are taken to 
be 3.4 A, the natural inter-layer spacing for graphite. 10 
In all cases, a semi-infinite vacuum layer exists above the 
top (z = 0) layer. 

Individual complex dielectric functions are obtained for 
each layer. We calculate the values through ab-initio 
methods for graphene and hBN. Empirical values from 
the literature are used for the metal and silicon substrate 
layers. These values are then stored for use as input to 
a continued-fraction algorithm, which yields the effective 
dielectric function. This algorithm is outlined in Sec- 
tion nm 



The effective dielectric function, £ (w, k, z), of the strat- 
ified structure in Fig. Q] is that of Lambin et al4 The 
expression for £ was derived from EELS theory in a re- 
flection geometry. The expression has been shown to be 
applicable to both phonons^ and polaritonsii in strat- 
ified structures with histogram-like dielectric functions 
(continuous within each layer) and interacting interfaces. 
Though Lambin's initial paper containing the derivation 
focused on semiconducting materials, the expression and 
the formalism from which it is derived also apply to sur- 
face plasmon structure of alternating of metal-insulator 
layers i 11 ! 12 It is worth noting that the specular reflection 
model agrees well with the spectroscopic predictions of 
the Bloch hydrodynamic model in the small wave vector 
regime considered in our work^ 

The z coordinate is in the direction perpendicular to 
the free surface of the sample, extending from the z = 
surface to — oo. k denotes the surface excitation (plas- 
mon or phonon) wave vector and ui is the frequency of 
excitation. 



£(k,ij,z) 



iD(k, 



E(k,w,z) ■ k/fc' 



(1) 



where D(k,cj,z) = e(to, z)E(k, uj, z), and e(ui,z) is the 
long wavelength dielectric function (tensor) of the mate- 
rial at z. £ remains continuous even in the case of sharp 
interfaces parallel to the x-y directions below the surface 
(as is the case in our multilayer system). This is due to 
the interface boundary conditions: continuity of D±_ and 

The effective dielectric function £$(k,u)) (Eq. [3]) is a 
solution to the Riccati equation (Eq. [5]), in the long- 
wavelength approximation k m 0, at the z — surface. 8 
We fix k as k = 0.05 A -1 for both the ab-initio calcula- 
tions and the composite calculation. Eq. [2] was derived 
for heterogeneous materials made of a succession of lay- 
ers (with homogeneous dielectric functions within each 
layer), the layers having parallel interfaces. e(z) are com- 
plex functions, with positive imaginary parts at z = 0. 8 



1 dg(z) 
k dz 



e(z) 



e(z) 
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and 



ai + a 2 — ja- 



ai = ej coth(fcc?i) 
hi = tij sinh(Mj). 
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Once individual dielectric functions are obtained, this 
procedure allows for the performance of mesoscale EELS 
calculations of a wide variety of structures. Layer thick- 
ness and materials can easily be substituted in the calcu- 
lation, with each EELS calculation running in a fraction 
of a second (nearly independent of the spectral range). 
EELS are calculated directly from the effective dielectric 
function as 



EELS = Im 



-1 



(6) 



Inspection of Eq.[3]reveals that for Im[eJ > 0, Im[£ ] > 
0. The EELS spectra given by Eq. [5] are then positive. 



C. Noble-Metal Dielectric Functions 

The copper, silver, and gold complex dielectric func- 
tions are empirical values by Johnson and Christy 1 ob- 
tained by reflection and transmission spectroscopy on 
vacuum-evaporated films at room temperature. Film- 
thickness in the Johnson and Christy study ranged from 
185 A- 500 A. Dielectric functions in the film-thickness 
range 250 A - 500 A did not vary significantly. In our 
work, the intermediate value 340 A was chosen to repre- 
sent bulk mode dominant (yet still nanoscale) metallic 
thin films. 



D. SiC>2 and Si Dielectric Constants 

Relative static permittivities of 3.9 and 11.68 are used 
for the Si0 2 and Si dielectric constants, respectively. 
These are reasonable and widely-used values from the 
literaturej 14 ' 15 



E. Graphene and hBN Individual Layer Dielectric 
Functions 



Complex dielectric functions for graphene and hBN 
are displayed in Fig. [5] (a) and Fig. [5] (b), respectively. 
These ab-initio calculations use the time-dependent den- 
sity functional theory in the local density approximation 
(LDA), and are implemented in the Python code GPAW, 
a real-space electronic structure code using the projec- 
tor augmented wave methodj 16 i 17 ~ — Both graphene and 
hBN dielectric functions are calculated in the armchair 
configuration for a momentum transfer value of 0.05 A -1 , 
along the T-M direction of the surface Brillouin zone. 
The armchair configuration and value of k are selected for 
comparison of EELS with existing data in the literature 
obtained via density functional theory (DFT) methods. 
The fc-point sampling with 30 x 30 x 1 Monkhorst-Pack 
grid was chosen for the band-structure and EELS calcu- 
lations for both graphene and hBN. 



Our model utilizes dielectric functions due to surface 
parallel excitations only, as the effective dielectric func- 
tion is derived in a specular reflection geometry. This is a 
reasonable approximation as out-of-plane excitations are 
minimal in graphene at energies less than rs 10 eV, and 
extreme UV radiation is outside of the regime of interest 
for this studyi^I A lattice constant of 2.46 A is used for 
graphene, and 2.50 A for hBN, as hBN is nearly isomor- 
phic to the graphene, except for the slightly larger lattice 
constant j2I 

The dielectric function we have obtained for graphene 
(see Fig. [2] (a)) is nearly identical to those of Yan et 
ali^ -, where the authors used the projector augmented 
wave methodology implemented in GPAW in terms of 
linear combinations of atomic orbitals, with a momen- 
tum transfer of 0.046 A -1 , along the T-M direction of 
the surface Brillouin zone. Both our results and that 
of Yan et al. display collective peaks for free-standing 
single-layer graphene at w 5 eV and w 15 eV. These val- 
ues also agree with experimental EELS results for single- 
layer graphene plasmons (with in-plane excitation), for 
example Eberlein et ali^2 find the tt plasmon at 4.7 eV 
and tt + a at 14.6 eV. 

A low-energy peak (below 1 eV) is also apparent in 
Fig. [2] (a). This feature, which has been observed in ab- 
intio calculations by other s 20 ' 24 , corresponds to a broad 
shoulder in the EELS, and is due to the low-energy tt 
— > tt* single-particle excitation. Low energy graphene 
features, including intraband transitions are generally 
quite dependent upon the value of momentum transfer^ 
These features are outside the energy regime of interest 
for the surface and bulk plasmon peaks in this paper, 
work in which very low-energy excitations are a focus 
should not neglect momentum transfer dependence. In 
the context of these calculations, increasing momentum 
transfer values lead to a shift towards higher energies of 
the tt plasmons for intrinsic graphene! 20 ' 26 Others, for ex- 
ample Gao et al. — have found (via the time-dependent 
local density approximation) linear dispersion for the tt 
plasmon (the key feature in our energy regime of inter- 
est) in single-layer graphene for both T-M and T-K di- 
rections. The calculations of Gao et al. were found to 
agree well with experiment. 

Our hBN dielectric function (see Fig. [21(b)) compares 
well with that of Yan et ali^i, who obtained their results 
in the long- wavelength limit, using the LDA-adiabatic 
local density approximation method. Both our results 
for hBN spectra and that of Yan et al. display a broad 
absorbtion peak with onset at « 4.5 eV and maximum 
at « 5.75 eV, which is in good agreement with previ- 
ous literature^ The dielectric function obtained for hBN 
is also similar to the experimental results of Tarrio and 
Schnatterly^ 9 -, where the authors measured a peak in the 
imaginary part of the dielectric function at w 6 eV. Dif- 
ferences between the complex dielectric function used in 
this study and that of 29 are attributable to the larger 
momentum transfer value of 0.13 A -1 in the Tarrio and 
Schnatterly study. 
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1 2 3 4 5 

energy (eV) 

FIG. 2. Complex relative dielectric function e(u>) for graphene 
(a) and hBN (b). Real and imaginary parts (e'(w) and e"(w)) 
are represented by solid and dotted lines, respectively. 



III. RESULTS 

A. Decreasing metallic film thickness and the 
begrenzung effect 

Figures [3] (a)-(c) (respectively, Ag, Au, and Cu on 
SiC^/Si, without graphene top coating) demonstrate the 
effect of decreasing noble metal layer thickness. These 
data should serve as a basis for comparison with Fig. 0] 
(a)-(c), which are discussed in detail in the next section. 

As the noble metal film thickness is reduced in Fig. [3] 
(a), the so called begrenzung effect is apparent. An in- 
crease in the surface-to- volume ratio in the metal causes 
enhanced coupling to the surface resonance and dimin- 
ished coupling to the bulk modesi 30 ' 31 

In the case of a thin metallic slab, empirical models 
have been quite thoroughly explored. Upon the introduc- 
tion of a boundary to an infinite metallic slab, a negative 
(begrenzung) peak is introduced at the same energy as 
the bulk peak, and a trailing surface peak appears, 3 ? The 
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FIG. 3. EELS: (a), (b), and (c) demonstrate the effect of 
differing thickness for the Ag, Au, and noble-metal layer, re- 
spectively, with a Si02/Si substrate, and without a graphene 
coating. Film thicknesses are 34 nm (solid line), 20 nm (long 
dashes), fOnm (short dashes), and 4nm (dotted line). 
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surface peak becomes more pronounced with decreasing 
thickness, as does the negative begrenzung peak, decreas- 
ing the net bulk-plasmon amplitude. A sharp transition 
can be observed, between 34 nm and 10 nm Ag film thick- 
ness, in which surface modes become dominant. This 
transition can also be observed for Au and Cu films, be- 
tween 34 nm and 20 nm. This is consistent with obser- 
vations in the well-validated and widely-used empirical 
data by Johnson and Christyji 



B. SPR Enhancement due to the graphene coating 

1. Si0 2 /Si Substrate 

Figures |31 (a) and (b) demonstrate the effect of the ad- 
dition of graphene coatings of varyious thicknesses to the 
Ag surface, in the case of a SiC^/Si substrate. For up to 
three layers of graphene, the surface peak remains virtu- 
ally unchanged. However, at three layers, the bulk peak 
amplitude is reduced by approximately 50%. This effect 
is attributed to the introduction of a thin boundary layer, 
decreasing the bulk resonance amplitude. Imposition of 
a boundary leads to simultaneous diminished coupling 
to bulk modes and enhanced coupling to surface modes, 
physically similar to the aforementioned begrenzung ef- 
fecti^Sr— At 20 layers, as we would expect, the bulk peak 
broadens significantly (indicating increased losses) and 
the system has spectral properties (broadening and peak 
position) resembling those of graphite on silver. Peak 
positions calculated for the 20-layer graphene coating 
(« 2.6 eV and s=s 3.5 eV) are very close to those measured 
in EELS of thin films of Ag nanoparticles evaporated 
on graphitic surfaces, for example 2.2 eV and 3.4eVi22 
Differences in peak position (particularly for the lower- 
energy peak) are attributable to differences in morphol- 
ogy between the Ag slab in our calculations and the film 
of Ag nanoparticles in the experimental samples. 



2. hBN Substrate 

Figure 0] (c) demonstrates the effect of the addition of 
graphene coatings of varying thickness to the silver sur- 
face, in the case of a hBN substrate. The use of the hBN 
substrate dramatically enhances the SPR peak, as well as 
shifting the surface peak to an energy very close to that 
of the bulk peak, even the case where no graphene coat- 
ing is used. This blue-shift of the SPR is due to the much 
smaller real part of the hBN dielectric function in com- 
parison to Si. Decreasing substrate dielectric function is 
known to dramatically blue-shift SPRi^i 

Addition of a single-layer graphene coating serves to 
further diminish the bulk peak without significant degra- 
dation of the surface peak. For up to 3 layers of graphene, 
the surface peak does not broaden significantly. As in the 
case of the Si02/Si substrate, broadening of the surface 
peak for more than 3 layers of graphene corresponds to 




3.0 3.5 

energy (eV) 




2.5 2.6 2.7 

energy (eV) 



c 

D 

>. 
i— 
ro 
i— 
+j 

!q 

03 
C 

o 

u 
c 

D 



o 




3.5 4.0 

energy (eV) 



FIG. 4. EELS for varying numbers of layers for the graphene 
film: (a) 34 nm Ag layer and S1O2/S1 substrate, (b) the sur- 
face peak in (a), and (c) 34 nm Ag layer and monolayer hBN 
substrate. In both cases the graphene layer numbers are 
(solid line), 1 (long dashes), 3 (intermediate-length dashes), 
10 (short dashes), and 20 (dotted line). 
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FIG. 5. EELS:(a) Comparison of the use of a Si02/Si sub- 
strate (solid line), a mono-, tri-, 10-, and 20-layer hBN sub- 
strate (long dashes, intermediate-length dashes, short dashes, 
and a dotted line, respectively). No graphene coating is used 
and the Ag film is 34 nm thick, (b) Close-up of EELS peaks 
for the hBN substrates. The suspended sample ( "vacuum sub- 
strate") is represented by the solid line in this case. Mono-, 
tri-, 10-, and 20-layer hBN substrates are represented by long 
dashes, intermediate-length dashes, short dashes, and a dot- 
ted line, respectively. 



increased losses. Indeed, at 10 graphene layers the re- 
sponse is what one would expect from a graphite coat- 
ing. A similar transition in plasmonic behavior at 10 
layers has been observed experimentally for multilayer 
graphene EELSi^ 



3. Ideal configurations for Ag: enhancement of surface 
modes 

Figure [5] depicts EELS for a 34 nm Ag film (with no 
graphene coating) on various substrates: SiC^/Si, mono- 



layer, tri-layer, and 10-layer hBN, and semi-infinite vac- 
uum (as would occur for a suspended sample). Introduc- 
tion of the high band-gap hBN substrate dramatically 
enhances surface modes as compared to SiC^/Si. The 
result is striking — the sample with a monolayer hBN 
substrate has a spectral profile nearly identical to that of 
a suspended sample. The EELS remain very nearly iden- 
tical to the suspended case for up to 10 layers. This is 
in accordance with experimental results that have found 
the electron mobility of graphene on hBN to be nearly 
that of suspended graphenei^ The similarity of the hBN 
substrate to the vacuum is expected due to its ultra- wide 
bandgap, and in the case of our model partly due to the 
atomic layer thickness. Upon inspection of Fig. [5] (b), 
it is clear that both real and imaginary parts of hBN's 
dielectric function are small; therefore at a thickness of 
3.4 A, the properties of this substrate approach those of 
the vacuum. The hBN substrate appears to dramatically 
enhance coupling to the surface plasmon resonance when 
compared to the SiC^/Si substrate, as well as reducing 
the peak strength of the lossy bulk plasmon. 

Fig. [5] (a) demonstrates the combined effects of 
graphene coating and hBN substrate for various Ag film 
thicknesses. Comparison of Fig. [6] (a) with Fig. [3] (a) re- 
veals that the introduction of a tri-layer graphene coat- 
ing and hBN substrate to the Ag film results in a surface 
plasmon intensity and relative surface to bulk intensity 
similar to a significantly thinner Ag film on the SiC^/Si 
substrate. For example, the surface plasmon intensity 
and relative surface to bulk intensity of the 34 nm Ag case 
in Fig. [5] (a) is comparable to the 20 nm case in Fig. [3] (a). 
Additionally, the surface plasmon intensity and relative 
surface to bulk intensity of the 20 nm Ag case in Fig. [5] 
(a) is comparable to the 10 nm case in Fig. 03a). 

Comparison of Fig. [6] (a) with Fig. [6] (b) (which shows 
spectra for free-standing Ag of various thicknesses) fur- 
ther elucidates the effect of a few graphene layers. For 
thicker Ag films, a few-layer graphene coating diminishes 
the bulk peak without significant broadening of the SPR. 
However, for very thin Ag films, of 4nm, for example, 
the effect of diminishing the bulk peak is negligible and 
the addition of the graphene coating of three layers only 
slightly broadens the SPR. This is in keeping with the 
idea that for bulk Ag, the addition of a graphene layer en- 
hances surface coupling, through the begrenzung effect. 
This is further demonstrated by the lack of significant 
broadening of the SPR peak for a single-layer graphene 
coating (see Fig. H(c)). 



C. Au and Cu Noble-Metal Film Composites 

1. Spectral changes of Au and Cu films due to the graphene 

coating 

One may wish to employ a noble metal other than 
Ag in multilayer structures such as depicted in this pa- 
per. Cu has the advantage of lower cost and is therefore 
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FIG. 6. (a)EELS: trilayer graphene on Ag layers (of vary- 
ing thickness) with a mono-layer hBN substrate, (b) Free- 
standing Ag films of various thicknesses. In both (a) 
and (b) Ag layer thicknesses are 50 nm (solid line), 34 nm 
(long dashes), 20 nm (intermediate- length dashes), and 10 nm 
(short dashes), and 4nm (dotted line). 

attractive for industrial applications. Both Cu and Au 
may also be of interest due to inherent surface plasmons 
that occur at longer wavelengths than those of Ag. In 
this section we demonstrate the effect of graphene coat- 
ings and the hBN substrate on Au and Cu films. For 
reference, Fig. [3] (b) and (c) demonstrate the effect of re- 
duction of noble-metal film thickness for Au and Cu films 
(respectively) on a SiC^/Si substrate. The begrenzung 
effect is again increasingly apparent in both Au and Cu 
film plasmonic response as film thickness decreases. At 
20 nm, both Au and Cu display strong SPR and broad 
absorption for higher energies, rather than a bulk plas- 
mon peak. 

Inspection of Fig. [7] (a) and (b) reveals that the intro- 
duction of a graphene coating to an Au or Cu surface 
does measurably strengthen surface modes while reduc- 



FIG. 7. EELS: (a) The effect of various numbers of layers 
for the graphene film, in the case of a 34 nm Au layer and 
SiC>2/Si substrate, and (b) The effect of varyious numbers 
of layers for the graphene film, in the case of a 34 nm Cu 
layer and SiC>2/Si substrate. Graphene layer numbers are 
(solid line), 1 (long dashes), 3 (intermediate- length dashes), 
fO (short dashes), and 20 (dotted line). 

ing bulk plasmon intensity, although the effect is not as 
dramatic as in the Ag case (see Fig. 0J. Coatings as 
thick as 20 layers further enhance and do not significantly 
broaden the surface peak. A jump in broad absorption 
(corresponding to increased losses) at higher energies is 
apparent in the case of both Au and Cu for more than 
three graphene layers. 



2. Metallic films with hBN substrate: noble-metal 
comparison 

Figure IS (a) and (b) are the EELS for Au and Cu 
films of thickness 34 nm (respectively), with monolayer 
graphene coating, comparing the Si02/Si substrate, hBN 
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FIG. 8. EELS: 34 nm Au (a) and Cu (b) films, with monolayer 
graphene coatings and various substrates including Si02/Si 
(solid line), hBN monolayer (dashes), and vacuum (dotted 
line) substrates. 



substrate, and suspended samples. Surface peaks at 
« 2eV for both Au and Cu are surface peaks due to 
change in dielectric function across the metal/SiO^/Si 
substrate interface. The addition of the hBN layer in 
these cases appears to sharpen and enhance Au and Cu 
inherent plasmons. Notably, in the Au and Cu cases, 
the hBN substrate appears to produce a similar spectral 
profile to the suspended "vacuum substrate", as in the 
case of the Ag film. Rigorous optimization of both the 
surface coating and substrate for both Au and Cu films 
is the subject of further investigation. 



IV. DISCUSSION 

In this study we systematically investigated the effect 
of the use of different plasmonic materials, different semi- 
conducting substrates, and different layer thicknesses in 



a multilayer graphene-based plasmonic composite struc- 
ture. For mono-, bi-, and tri-layer graphene, in the case 
of 34 nm Ag layers (where bulk modes would dominate 
in the absence of a graphene coating), bulk plasmon 
modes are significantly diminished, while maintaining the 
strength of surface modes. 

Reduced plasmon losses for graphene coated plasmonic 
metals on hBN substrates likely originates within the 
context of our mesoscopic model from two primary physi- 
cal effects: (1) The addition of a graphene boundary layer 
on the metallic surface reduces coupling of excitations to 
bulk plasmons through the begrenzung effect. The origin 
of the begrenzung effect is a reduction of the degrees of 
freedom for excitations, and thus further surface confine- 
ment comes at the expense of bulk oscillations, leading 
to reduced losses. (2) The strongly insulating hBN sub- 
strate diminishes bulk losses and enhances surface con- 
finement particularly through a reduction in scattering. 

Traditional substrates such as Si are known to degrade 
the electron transport properties of graphene as com- 
pared to suspended samples. - Reduced transport prop- 
erties are due to various scattering mechanisms. There 
is strong evidence that scattering in these systems is due 
in large part to various substrate interactions including 
intcrfacial phonons (which are taken into account in our 
work), surface charge traps, and substrate stabilized rip- 
ples.— 

It is therefore reasonable to mitigate scattering in 
graphene-based systems for plasmonics (and devices for 
plasmonics in general) by suspending samples, or through 
the use of alternative substrates^ In this work, for Ag 
films, the use of a hBN substrate (rather than the more 
traditional silicon) is found to shift plasmonic coupling 
towards surface modes, both diminishing bulk losses and 
enhancing the SPR peak amplitude. This effect is par- 
ticularly dramatic in the case of the 34 nm Ag film (in- 
cluding those with no graphene top-coating), as an Ag 
film at this thickness on an Si02/Si substrate produces 
EELS that are extremely dominated by the bulk reso- 
nancei This is interesting in light of recent experimental 
results, showing electron mobility of graphene on hBN 
to be similar to that of suspended graphene^ We expect 
SPR enhancement for other wide-bandgap substrates due 
to enhanced surface plasmon field confinement. Further 
comparison of alternate wide-bandgap substrate materi- 
als, such as SiC, are the subject of future work. 

As expected, bulk modes were quenched for 4nm Au 
and Cu films. For Ag, bulk modes very nearly vanished 
at this thickness, and are nearly undetectable when cou- 
pled with the hBN substrate. However, the morphol- 
ogy of very thin metallic films deposited on graphene is 
difficult to control, often forming island-like structures 
of various sizes 4 This is also the case for graphene de- 
posited on a metallic substrate, as agglomeration below 
a critical thickness is a general property of thin films. 4 
Results of this study indicate that Ag films as thick as 
34 nm, when coated with a graphene film and placed on 
a wide-bandgap substrate such as hBN, may also be em- 
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ployed for low-loss plasmon resonance applications. An 
hBN substrate, due to its ultra-flat morphology, would 
also lend additional uniformity to the structure. 

The mesoscopic model used in these calculations has 
several limitations that are worth discussing. Results of 
this study are valid in the long- wavelength limit for which 
the continued fraction expression by Lambin et al. was 
derived. Additionally, coupling between layers is clas- 
sical (via boundary conditions), and as a result inter- 
layer hopping is neglected. This tunnelling has been de- 
termined by angle-resolved photoemission spectroscopy 
to be important for graphene band-structure, leading 
to 7r band splitting, which increases with layer number 
(w 0.7 eV at 4 layers)^ Multilayer graphene excitation 
spectra have been calculated by Ohta et al. with the in- 
clusion of inter-layer hopping in the kinetic Hamiltonian. 
Tunneling is found to be primarily important in accu- 
rately producing the low energy region of the excitation 
spectrum (lu ~ Vpk), where vp is the Fermi velocity) i2£ 
Though tunnelling is not as critical in the regime of inter- 
est for this work, calculations for graphene-based multi- 
layer plasmonic systems geared towards lower-energy ap- 
plications should take care to include inter-layer hopping. 

More general EELS, accurate for a wider range of k 
values and angles of incidence will require first princi- 
ples calculations for the entire structure, and incorporate 
higher level quasiparticle interactions (including the ef- 
fect of excitons). Using guidance from the results of this 
study, future work consisting of 1-3 layers of graphene on 
silver, with the hBN substrate are the subject of future 
work. The band-gap for hBN is underestimated by about 
33% in the LDA22 — the effect of hBN on surface con- 



finement of plasmons may be even more dramatic than is 
demonstrated in this work. A quasiparticle GW correc- 
tion to the LDA calculation will be employed, as it has 
been demonstrated to bring the hBN bandgap into close 
agreement with experimental results. 39 

As the individual layers in these structures are nearly 
isomorphic, rather than exactly isomorphic, future work 
will incorporate the effect of lattice strain on the opti- 
cal properties of the composite. Strain engineering is 
expected to provide a further means of plasmon res- 
onance tunability42 Additionally the effect of config- 
urations other than the flat-armchair morphology for 
graphene is the subject of ongoing work. The effect of 
rippling and defects in graphene-based multilayer struc- 
tures is of particular interest, as is a rigorous quantifica- 
tion of the uncertainties involved in the modelling and 
manufacturing stages, and optimal design against those 
uncertainties. 
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